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We study the electromagnetic fl c 7 —> fi* transition in 2+1 flavor lattice QCD, which gives access 
to the dominant decay mode of fl* baryon. The magnetic dipole and the electric quadrupole 
transition form factors are computed. The magnetic dipole form factor is found to be mainly 
determined by the strange quark and the electric quadrupole form factor to be negligibly small, in 
consistency with the quark model. We also evaluate the helicity amplitudes and the decay rate. 

PACS numbers: 14.20.Lq, 12.38.Gc, 13.40.Gp 


I. INTRODUCTION 

Recently there has been a significant progress in our understanding of the heavy-flavor hadron sector. Experimen¬ 
tally, all the ground-state single-charmed baryons and several excited states, as predicted by the quark model, have 
been confirmed pQ. Unlike other single-charmed baryons, a precise observation of fl c baryon had been long overdue. 
Only recently, Belle Collaboration has made a rigorous experimental study of f l c using the decay fl c ° -+ Q~n + 3]. 

The 12/(css) has the quantum numbers J p = | + and is the heaviest known single charmed hadron that decays 
weakly. Within the multiplet structure of flavor SU(4), Q c belongs to a sextet of flavor symmetric states, which sits 
on the second layer of the flavor mixed-symmetric 20-plet. The average mass value reported by the Particle Data 
Group (PDG) is m Qc = 2695.2 ± 1.7 MeV [I]. 

The excited Q*°(css) baryon was first observed by BABAR Collaboration in the radiative decay Q* —> fl c 7 [3], 
Belle Collaboration has confirmed their observation by reconstructing Ul* in the same radiative decay mode m- Th ey 
measured the relative mass difference with respect to the ground state mn* — mn c = 70.7 ± 0.9+Q.g MeV in very good 
agreement with the BABAR observation. The average mass value reported by PDG is inn* = 2765.9 ± 2.0 MeV pQ- 

The quantum numbers have not been measured but natural assignment is that it completes the ground state J p = | + 
sextet, which sits on the second layer of the flavor-symmetric 20-plet of SU(4). The mass difference with respect to 
the ground state is too small for any strong decay to occur, therefore the radiative channel f l* —> Ul c 7 is the dominant 
decay mode. 

In addition to BABAR and Belle Collaborations, experimental facilities such as LHCb, PANDA, Belle II, BESIII 
and J-PARC are expected to give a more precise determination of charmed baryon spectroscopy. Concurrently, recent 
lattice-QCD studies on the spectroscopy of charmed hadrons are also very promising. The ground-state charmed 
baryons with spin up to 3/2 have been studied in quenched |U[S] and full QCD [BHS]. The results for baryon masses 
as determined from lattice QCD are in good agreement with experiment. 

Recently we have examined the charmed baryons in lattice QCD in order to reveal their electromagnetic struc¬ 
ture mm- We have extracted the charge radii and magnetic moments of J = 1/2 charmed baryons by computing 
their elastic electromagnetic form factors on the lattice. A similar study for J = 3/2 baryons is in progress. A phe¬ 
nomenologically more interesting problem is the electromagnetic transitions between J = 1/2 and J = 3/2 baryons, 
which are more accessible by experiments as explained above. 

Being motivated by the experimental discovery of the f2* baryon in the radiative decay mode, in the present work we 
give a timely study of the J = 1/2 —>■ J = 3/2 electromagnetic transition of single charmed strange baryons in lattice 
QCD. In particular we study the fU 7 -+ fl* transition, which gives access to three Sachs form factors, the helicity 
amplitudes, the decay width and the lifetime. This work is reminiscent of Refs. fUjHTo] . where the electromagnetic N to 
A, and the other octet to decuplet transitions have been studied. The electromagnetic transitions of charmed baryons 
have also been studied within heavy hadron chiral perturbation theory |16H18j and in the quark models nuu. 

The three transition form factors, namely, the magnetic dipole (Ml), the electric quadrupole (E 2) and the electric 
charge quadrupole [C 2) provide valuable information about the structure and shape of J = 1/2 and J = 3/2 baryons. 
Earlier studies have focused on the transition moments between N and A. Experimentally, pure single spin-flip 
Ml transition has been found to dominate. Of special interest is the small but non-vanishing values of E2 and C 2 
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moments, implying that the shapes of N and A deviate from spherical symmetry f22]. Quark models predict a nonzero 
value for E2 and C 2 [23], which has also been confirmed experimentally [24] [25]. However, the results from various 
theoretical approaches are not in complete agreement with experiment and this issue is still unsettled. 

The experimental results for fl c 7 —> f l*, on the other hand, are not yet precise enough to allow a determination of 
the transition strengths. In this work, we will mainly focus on the Ml and E 2 transition form factors. Unlike in the 
case of Nj —> A, the mass splitting between fl c and fi* can be reproduced on the lattice and an accurate contact can 
be made to phenomenological observables via these two form factors. We employ near physical 2+1-flavor lattices that 
correspond to a pion mass of approximately 156 MeV. The data for electromagnetic transition form factors are often 
noisier than those for elastic form factors, particularly for C 2 form factor. Considering also the limited number of 
gauge configurations we have at the smallest quark mass, we study the Ml and E 2 form factors for the lowest allowed 
lattice momentum transfer. We, however, make contact with the transition moments at zero-momentum transfer by 
assuming a simple scaling at low momentum transfer values [12] . 


II. LATTICE FORMULATION 

Electromagnetic transition form factors for O c 7 —> O* can be calculated by considering the baryon matrix elements 
of the electromagnetic vector current = |c(x)7 /1 c(x) — ^s(x)j /1 s(x). The matrix element can be written in the 
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following form 


(Q:(pV)|^|O c (;P, S )) 


. j2 f TO* TO ^ 

*V3 U*(P')£(P)/ 


Ur(p\s')0 Tll u(p,s), 


with the operator O tM given in terms of Sachs form factors as [26] 

O t » = G M1 (q 2 )Kl £ + G E2 {q 2 )KZ + G C 2 (q 2 )K T c », 

where 


(1) 

( 2 ) 


K lii = -3^(m* + to) 2 - q 2 ^j ie T>Jl (Pq) (to* + to)/2to, (3) 

K E 2 = - 60 _1 (g 2 ) ie Tp {Pq) e^(p'q) 7 5 (to* + to)/to, (4) 

K Ch = -30 _1 (g 2 ) q T {q 2 P> 1 - q ■ P q 11 ) «7 5 (to* + m)/m. (5) 

Here p and p' denote the incoming and the outgoing momenta, respectively, q = p 1 — p is the transferred four- 
momentum, P = (p' +p)/2 and 


U(g 2 ) = (( TO * + m ) 2 - 9 2 ) ((m* - to) 2 - q 2 ^ . (6) 

We use the shorthand notation e TIX (Pq) = e T/ ‘ QI ' P a q u . The spins are denoted by s, s' and the masses of U* and fl c by 
to* and to, respectively. u(p , s) is the Dirac spinor and u T (p , s ) is the Rarita-Schwinger spin vector. For real photons, 
G c 2(0) does not play any role as it is proportional to the longitudinal helicity amplitude. 

The Rarita-Schwinger spin sum for the spin-3/2 held in Euclidean space is given by 


s)u T (p, s) 


—i'y ■ p + m* 
2 TO* 


1 

9 <tt ~ 2 + 


2 PaPr 
3 TO* 


Pair ~ Prig 
3 TO* 


( 7 ) 


and the Dirac spinor spin sum by 


Y,u(p,s)u(p, s ) 

S 


— i'y ■ p + to 
2 TO 


( 8 ) 


We refer the form factors Gm 1, G E 2 and Gc 2 as the magnetic dipole, the electric quadrupole and the electric charge 
quadrupole transition form factors, respectively. 

To extract the form factors we consider the following matrix elements, 

(GaT Qc (I; p; T 4 )) = 5>- p ' x rr' x (vac|T[77“(z)?7“'(0)]|vac), 


X 


( 9 ) 
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<G nA (i;p;r 4 )> = 5>" <P ' Xr 4 a ' x <vac|T[?f (aOif'(0)]|vac), 

X 


{G a oJ c (t 2 ,t 1 ;p\p;T)) =-i 


3 -'P'X2„iq'Xii 


( v ac|r[7?“(®2)j M (®i)?? Q (0)]|vac), 


with the spin projection matrices defined as 


r < = 


1 

2 



r 4 


1 

2 



( 10 ) 

( 11 ) 


( 12 ) 


Here, a, a' are the Dirac indices, a and r are the Lorentz indices of the spin-3/2 interpolating field and cr, are the 
Pauli spin matrices. An initial f2 c state is created at time zero and interacts with the external electromagnetic field 
at time t\. At time t 2 the final fl* state is annihilated. 

The baryon interpolating fields are chosen, similarly to those of A and N as 


= -^f 3k {‘2[s T \x)C'y l _ l c j (x)\s k (x) + [s T \x)C-/„s J {x)]c k {x)} , 


r)(x) = e l:>k [s Tl (x)Cj 5 c , (x)\s k (x), 


(13) 

(14) 


where i, j, k denote the color indices and C = y 4 72- It has been shown in Ref. J5] that the interpolating field in 
Eq. (13) has minimal overlap with spin-1/2 states and therefore does not need any spin-3/2 projection. 

To extract the form factors, we calculate the following ratio of the two- and three-point functions: 


p',p;r ; /x) = 


<g; 


o*j»‘n c 


(t 2 ,ii;p',p;r)) 


(^G^(f 2 ;p';r 4 )> 


5 w Gg :n '(2t 1 ;p / ;r 4 ))l V2 


QSi c n c ( 2 f 4 ; p; r 4 )) 


(15) 


In the large Euclidean time limit, t 2 — t\ and t\ a, the unknown normalization factors cancel and the time 
dependence of the correlators can be eliminated. Then the ratio in Eq. (15) reduces to the desired form 


Rcr(t2,ti-, p',p;T;n) 


ti»q 


> n CT (p',p;T;/Li), 


(16) 


where we can look for a plateau to extract the form factors. We choose the ratio in Eq. (15) among several other 
alternatives used in the literature OGHISj due to the good plateau region and the quality of the signal it yields. 

We single out the Sachs form factors by choosing appropriate combinations of Lorentz direction /r and projection 
matrices T. When S2 C is produced at rest and momentum is inserted in one spatial direction, we have |13| 


G C 2 (q 2 ) = G(q 2 )^n fc (q, 0; iT k - 4) 


Gmi(<? 2 ) = C(q 2 ) 1 


G E 2 (q 2 ) = C( q 2 ) 1 


m* 


nz(<?fc,0;r fe ;Z) - — U k (q k , 0; T;; l) 

-C/* 


TO* 


^i(q k , 0; T fe ; l) + — n fc (g fc , 0; T ( ; l) 

-C/* 


where 


C(q 2 ) = 2v / 6^^- 

to + TO* 


1 + 


TO*y /2 / + ^_y /2 

Ej \ + 3 ml) ' 


(17) 

(18) 

(19) 

( 20 ) 


and k and l are two distinct indices running from 1 to 3. When f2* is produced at rest, to* = EC* in Eqs. ( 17]|T9 | and 

2 \ 1/2 


G(q 2 ) = 2^6- 


Em 


( m\ 1 / 2 ( 

( 1+ E ) ( 1 + 


3 to 2 


( 21 ) 


We have run our simulations on gauge configurations generated by PACS-CS collaboration m with the nonpertur- 
batively 0(a)-improved Wilson quark action and the Iwasaki gauge action. The details of the gauge configurations 
are given in Table |l| The simulations are carried out with near physical u,d sea quarks of hopping parameter n u,d = 
0.13781. This corresponds to a pion mass of approximately 156 MeV m- The hopping parameter for the sea s quark 
is fixed to Kg ea = 0.13640 and the hopping parameter for the valence s-quark is taken to be the same. 
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TABLE I. The details of the gauge configurations used in this work m- We list the number of flavors ( Nf ), the lattice 
spacing (a), the lattice size ( L ), inverse gauge coupling (/3), clover coefficient ( csw ), number of gauge configurations employed 
and the corresponding pion mass 


w 3 X N t 

Nf 

a [fm] 

L [fm] 

P 

Csw 

# of conf. 

[MeV] 

32 3 x 64 

2+1 

0.0907(13) 

2.90 

1.90 

1.715 

194 

156(7) (2) 


As we perform the simulations at only one (near-physical) quark mass, a chiral extrapolation cannot be made. 
However, we can make an estimation of an uncertainty anticipated from such an extrapolation, based on our simu¬ 
lations of elastic f l c electromagnetic form factors. We have performed the chiral extrapolations for electric/magnetic 
charge radii and the magnetic moment of baryon in Ref. nn again, including the data at m^ ~ 156 MeV. We tried 
constant, linear and quadratic fit functions. For all cases, the chiral-extrapolated values and those at the smallest 
pion mass are in very good agreement within their error bars. Different fit forms we use imply a systematic error of 
less than 1%. Hence, we anticipate to have a similarly negligible error from such an extrapolation of Ml and E2 form 
factors here. 

For the charm quarks, we apply the Fermilab method [28] in the form employed by the Fermilab Lattice and MILC 
Collaborations [2SISD|. A similar approach has been recently used to study charmonium, heavy-light meson resonances 
and their scattering with pion and kaon [5TH33| . In this simplest form of the Fermilab method, the Clover coefficients 
ce and eg in the action are set to the tadpole-improved value 1 /mq, where uq is the average link. Following the 
approach in Ref. [3lj, we estimate uq to be the fourth root of the average plaquette. We determine the charm-quark 
hopping parameter k c nonperturbatively. To this end, we measure the spin-averaged static masses of charmonium 
and heavy-light mesons and tune their values accordingly to the experimental results, which yields k c = 0.1246 m- 

We make our simulations with the lowest allowed lattice momentum transfer q = 2ir/N s a, where N s is the spatial 
dimension of the lattice and a is the lattice spacing. This corresponds to three-momentum squared value of q 2 = 
0.183 GeV 2 . In order to access the values of the form factors at Q 2 = — q 2 = 0, we will apply the procedure in Ref. PJj 
and assume that the momentum-transfer dependence of the transition form factors is the same as the momentum 
dependence of the H* baryon charge form factor. Such a scaling is also suggested by the experimentally measured 
proton form factors and it was used in previous analyses such as baryon octet to decuplet electromagnetic transition 
form factors m- While extrapolations in finite momentum suffer from large statistical errors since one has to rely 
on a functional form, the scaling approach provides a more precise determination of the form-factor values at zero 
momentum transfer. In applying this procedure, we consider s and c quark sectors separately as their contributions 
to the charge form factors scale differently. For instance, the scaling of Gmi is given by 


G S m C i(0)=G^( 9 2 


Gg C (0) 

g s e c (q 2 Y 


( 22 ) 


The heavy-quark contribution yields a harder form factor whereas the light-quark contribution is soft and the form 
factor falls off more rapidly muon- The form factors are extracted in two kinematically different cases. In the first 
case, the f 1 * is produced at rest and the has momentum — q and in the second case, the H c is at rest and carries 
momentum q. 

In order to increase statistics, we insert positive and negative momentum in one of the spatial directions and make 
a simultaneous fit over all available data. We also consider current along all spatial directions. The source-sink time 
separation is fixed to 1.09 fm = 12a), which has been shown to be sufficient to avoid excited state contaminations 
for electromagnetic form factors HU. Using translational symmetry, we have employed multiple source-sink pairs by 
shifting them 12 lattice units in the temporal direction. All statistical errors are estimated by the single-elimination 
jackknife analysis. We consider point-split lattice vector current 

= 1/2 [q(x + n)U^(l + 7 ^)g(x) - q{x)U^(l - 7 ,,)q(x + n)], (23) 

which is conserved by Wilson fermions. 

A wall-source/sink method [34] has been employed, which provides a simultaneous extraction of all spin, momentum 
and projection components of the correlators. The gauge non-invariant wall source/sink requires fixing the gauge. We 
fix the gauge to Coulomb, which gives a somewhat better coupling to the ground state. The delta function operator 
is smeared over the three spatial dimensions of the time slice where the source is located, in a gauge-invariant manner 
using a Gaussian form. In the case of s quark, we choose the smearing parameters so as to give a root-mean-square 
radius of (ri) ~ 0.5 fm. As for the charm quark, we adjust the smearing parameters to obtain (r c ) = (ri)/ 3. 
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TABLE II. The fl c and fl* masses (at a pion mass of = 156 MeV) together with the experimental values [T] and those 
obtained by PACS-CS [6] (at the physical point). We have also included results by ETMC [9j and Briceno et al. [8] obtained 
by chiral extrapolation. 



This work 

PACS-CS [S] 

ETMC [U 

Briceno et al. [H] 

Exp. T| 

m [GeV] 

2.750(15) 

2.673(17) 

2.629(22) 

2.681(48) 

2.695(2) 

m* [GeV] 

2.828(15) 

2.738(17) 

2.709(26) 

2.764(49) 

2.766(2) 


III. NUMERICAL RESULTS AND DISCUSSION 


We extract the U c and 0* masses using the two-point correlators in Eqs. and ( [To] ) . Our results for the U c and O* 
masses are given in Tabic |U] together with the experimental values and those obtained by other lattice collaborations. 
While we see a few percent discrepancy between our results obtained at a pion mass of m n = 156 MeV and those of 
PACS-CS obtained at the physical point, the mass splitting to* — to is nicely produced in agreement with experiment. 
We define the sum of all correlation-function ratios as 


|q| 6 ^ 


n 2 = 


C(q 2 


x y^n fc (g fc ,0;r;;l), 


k,l 


so that Eq. (18) and Eq. (19) becomes, 


G M1 ( q 2 ) = n 1 - r ^n 2 , 
G E 2(q 2 ) = n 1 + ^n 2 . 


(24) 

(25) 

(26) 


Fig- [T| illustrates the III and n 2 as functions of the current insertion time, p, for s- and c-quark sectors separately. 
The two ratios have opposite sign and they add constructively when they are subtracted. We extract the form factors 
by fitting the correlation-function ratios by a horizontal line where a plateau develops. We illustrate both kinematical 
cases giving consistent results within their error bars. A clear plateau can be realized in both kinematical cases, being 
more flat when f l c is produced at rest. We fit the correlation function ratios in the range t\ = [3,6]. The statistical 
errors, on the other hand, are smaller when ff* is at rest. The values of the form factors from the two kinematical 
cases are consistent with each other. 

It is straightforward to extract G E 2 once we construct the correlation function ratios for Gmi- The correlation 
functions have opposite signs and are of similar magnitudes, which result in a vanishing value for G E 2 when they are 
added. We determine G E 2 by fitting Ip and n 2 separately and combining the results. This procedure gives consistent 
results with fitting the sum of the correlation ratios. 


Our numerical results are reported in Table III We give the values of Gmi and G E 2 form factors at the lowest 
allowed momentum transfer and at zero momentum transfer for the two kinematical cases as explained above. The 
quark sector contributions to each form factor are given separately. The form factors can be inferred from individual 
quark contributions by 


Gmi(Q 2 ) = l x G C M1 (Q 2 ) - i x G S M1 (Q 2 ), 


(27) 


and similarly for G^ 2 (Q 2 ). The values of the form factors at Q 2 = 0 are extracted using the scaling assumption in 
Eq. •|22li. 


Similarly to what has been observed in the case of elastic form factors m , M 1 form factor is dominantly determined 
by the contribution of the s-quark sector, which is approximately one order of magnitude larger than that of the c- 
quark sector. This pattern is consistent with hyperon transition form factors m The heavier quark contribution is 
systematically smaller than that of the light quarks. From a quark-model point of view, the coupling of the photon to 
the light quarks prevails in the heavy-quark limit and the heavy quark acts as a spectator. In this limit, the transition 
proceeds dominantly through the spin flip of the light degrees of freedom and only Ml transition is allowed. Only 
finite mass effects of the heavy quark may lead to a nonzero value of E 2 form factor. Our results show that the 
two quark sectors contribute with opposite signs and yield a value with a statistical error of approximately 5% when 
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rest. 


The correlation function ratios IIi and II 2 in Eq. (241 as functions of the current insertion time (ti) for s- and 
We also display GJj/ obtained using Eq. (251. The squares (triangles) denote the kinematical case when 12* 


c-quark 
(fi c ) at 


TABLE III. The results for Gm\ and Ge 2 form factors at the lowest allowed four-momentum transfer and at zero momentum 
transfer for the two kinematical cases. The quark sector contributions to each form factor are given separately. Note that the 
statistical uncertainty is large in Ge 2 and results are consistent with zero. 


Q 2 [GeV 2 ] 

G|n(Q 2 ) 

g c m1 (Q 2 ) 

Gmi(Q 2 ) 

g s E2 (Q 2 ) 

G c E2 (Q 2 ) 

Ge2(Q 2 ) 

12* at rest °' 18 ° 

1.257(67) 

-0.167(33) 

-0.530(28) 

0.041(132) 

0.008(26) 

-0.008(50) 

0 

1.622(87) 

-0.175(34) 

-0.657(33) 

0.052(171) 

0.009(27) 

-0.012(62) 

12 c at rest 0.168 

1.269(177) 

-0.174(37) 

-0.539(78) 

-0.035(124) 

0.061(25) 

0.052(48) 

0 

1.637(229) 

-0.183(39) 

-0.667(96) 

-0.045(160) 

0.064(27) 

0.058(60) 


combined via Eq. (271. The values from the two kinematical cases are consistent with each other within their error 
bars. 


In contrast, the extracted values of Gei at finite and zero momentum transfer are small and consistent with zero 
within their error bars. A comparison of Gmi and Ge 2 reveals that the transition is entirely determined by Ml 
transition. In quark model, the quadrupole transition moments arise from the tensor-induced D-state admixtures of 
the single-quark wavefunctions |23j and the two-quark exchange currents [351I3S]- In the first, the spins of the quarks 
remain the same but an S-state quark is changed into a D-state. The latter can be interpreted as the spin flip of a 
diquark inside the baryon. Given the dependence of the tensor force on the inverse quark mass, one would expect 
to obtain a smaller Ge 2 value for heavy baryons as compared to that in the light-baryon sector, in consistency with 
what we have found. The smallness of the E 2 form factor can also be understood as a chiral suppression. The E 2 
amplitude is dominated by pion loops and the leading contribution comes from chiral logs which can be computed in 
heavy-baryon chiral perturbation theory |37l l38| . 

The Sachs form factors calculated above can be related to phenomenological observables such as the helicity am¬ 
plitudes and the decay width. The relation between the Sachs form factors extracted in this work and the standard 
definitions of electromagnetic transition amplitudes /mi and /e 2 in the rest frame of D* are given by E1I13 


fMi(q 2 ) = 

\J4-ku I 


\ 1/2 G M i(q 2 ) 

2 TO ' 

^ m 2 

) [1 — q 2 /{m + m*) 2 ] 1 / 2 ’ 

fE2(q 2 ) = 

\J4-Ka I 

f \q\m * N 

\^ 2 GE2(q 2 ) 

2 m ' 

v rn , 

1 [l-g 2 /(m + m*) 2 ] 1/2 ’ 


(28) 

(29) 


where a = 1/137 is the fine structure constant. The helicity amplitudes A\i 2 and A 3 / 2 can be deduced from the 
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TABLE IV. The results for the helicity amplitudes and the decay width in the rest frame of fl*. The helicity amplitudes are 
given at finite and zero momentum transfer. The zero-momentum values are obtained using the scaling assumption in Eq. (221. 


Q 2 

/mi 

fE2 

Al/2 

A3/2 

r 

[GeV 2 ] 

10 -2 [GeV -1/2 ] 

10~ 2 [GeV 1/2 ] 

10~ 2 [GeV 1/2 ] 

10~ 2 [GeV 1/2 ] 

[keV] 

0.180 

0 

-0.795(42) 

-0.988(50) 

-0.012(75) 

-0.018(93) 

0.416(116) 

0.521(145) 

0.678(71) 

0.840(88) 

0.074(8) 


transition amplitudes as follows: 


^ 1/2 (? 2 ) — — 1/2[/mi(? 2 ) + 3/s2(<? 2 )]) 

^3/2(<7 2 ) = -\/3/2[/ mi (<7 2 ) - fE2 (*? 2 )]■ 


Then the decay width is given by [T] 


(30) 

(31) 


2 \ 2 

Y_m*mi m \ nA 


8n 


{l^i/2 (0)p + |A 3/2 (0)f }, 


(32) 


where we have used the constraint q = (m 2 — ?n 2 )/2m* at q 2 = 0. The decay width can also be obtained from the 
Sachs form factors: 


r=S^ !! ^F{3|c ra (0)P + | GM1 (0)iq. 


(33) 


Since the above formulas are continuum relations, we use the experimental values of Q c and masses in calculating 
the helicity amplitudes and the decay width. Our numerical results for the helicity amplitudes in the rest frame of 
O* and the decay width, at finite and zero momentum transfer, are reported in Table IV A comparison to the 
N 7 —► A transition (T| reveals that, the helicity amplitudes are suppressed roughly by five orders of magnitude due 
to diminishing contribution of the heavy quark, the overall reduction in the transition form factors and the larger 
baryon masses. 

Since no strong decay channel is kinematically allowed, the total decay rate of Q* is almost entirely in terms 
of the photon decay mode. Eventually a significantly suppressed value of the fi*-baryon decay width is yielded, 
making Q* one of the longest living spin-3/2 charmed hadrons. The suppression in the decay width can be mainly 
attributed to the small f2*-f2 C mass splitting. The decay width in Table IV is translated into a lifetime of r = 1/T = 
8.901(913) x 1CT 18 sec. 

The electromagnetic transitions of charmed baryons have also been studied within heavy hadron chiral perturbation 
theory PHUT] and quark models [T9l - i2l] . It has been found that the charmed baryon electromagnetic decays are 
suppressed, in qualitative agreement with our result. Of special interest is the E* ,+ —> E+q decay having a similarly 
small width in the quark model m- An enhanced width is foreseen in the E* + —» A+q decay, which would be 
interesting to study on the lattice. The literature on n c q —> D* transition is limited. Non-relativistic quark model 
prediction for S2* decay width [19| is one order of magnitude larger than what we have calculated in this work. Note 
that given the small L 2 *-D c mass splitting, such a large width would require a Gmi value as large as that of Ay — » A 
transition. This cannot be justified as we have found that the heavy-quark contribution diminishes and there is no 
indication that the light quark contribution is enhanced. 

In conclusion, we have computed the f2 c q —» Q* transition in lattice QCD. The dominant contribution is due to 
the magnetic dipole form factor, which we have calculated with a statistical precision of about 5%. The electric 
quadrupole transition has been found to be negligibly small in consistency with the quark model. We have extracted 
the helicity amplitudes and the decay width, which have been found to be suppressed. This transition is of particular 
interest because of its relevance to current and proposed experimental facilities such as LHCb, PANDA, Belle II, 
BESIII and J-PARC, which are expected to measure the electromagnetic decay widths of charmed baryons with a 
higher precision. 
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